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Abstract 

One of standard methods to predict the phenomena of squeezing consists in splitting the 
unitary evolution operator into the product of simpler operations (Yuen [1], Ma and Rhodes 
[2]). The technique, while mathematically general, is not so simple in applications and leaves 
some pragmatic problems open. We report an extended class of exponential formulae, which 
yield a quicker insight into the laboratory details for a class of squeezing operations, and 
moreover, can be alternatively used to programme different type of operations, as: 1) the 
free evolution inversion, 2) the soft simulations of the sharp kicks (so that all abstract results 
involving the kicks of the oscillator potential, become realistic laboratory prescriptions). 


1 The manipulation problem 

Below, we shall dissent from the orthodox subject of “squeezed states” and dedicate some attention 
to a more general problem. Suppose, one has a quantum system whose states are represented by 
vectors in a Hilbert space H. Now, choose any unitary operator 

U :7i ->H (1) 

Can U be achieved as a realistic evolution operation, performed under the influence of some 
external fields? 

The problem so stated, belongs to the quantum manipulation theory , a domain which has 
progressed quickly in the last decades. The first cases of the dynamical manipulation (for a finite 
dimensional space of states) achieved wide publicity under the name of the spin echo (e.g. [3]). 
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The general problem of manipulation (control) of quantum states dates from the works of Lamb 
[4], Lubkin [5] and followers [6, 7, 8. 9, 10). Quite independently, the subject has been put forward 
in quantum chemistry where it may soon become crucial (11, 12). For an infinite dimensional 
7 ■{ = L 2 (R) some dynamical operations present a considerable challenge but only one of them has 
become a “conference subject” . We of course refer to the operation of squeezing-. 

U = e (* at2-z * ai )/ 2 (general squeezing) z G C (2) 

and/or 

U = e lX ( qp+pq ^' 2 (scale transformation, A G R (3) 

coordinate squeezing) 

Note, that there are several concepts of squeezing in the literature. By choosing (2-3) we ask 
about the “operatorial squeezing”, i.e. the shape transformation which affects all wave packets 
alike, independently on their initial form. Thus, under the influence of (3) the canonical observ- 
ables q,p are transformed into 

U\U = e~ x q 

U^pU = e x p (4) 

and simultaneously all the wave packets ip = {i>(x)} are deformed as: 

(Uip){x) = Vk ip(kx), k = e x (5) 

As found by Yuen [1], the simplest method of producing such effects in L 2 ( R) consists in 
application of variable oscillator potentials with the time dependent Hamiltonians: 

tf(0 = y + ^(0 2 |- ; [?.pI = * ( 6 ) 

and the most explicit illustrations of this fact can be found in the exponential formulae, which 
express the evolution operator U(t) (generated by (6)] as the product of simpler exponential 
operations. 

The very subject of the exponential identities has already some antiquity, starting from the 
papers of Zassenhaus, Baker, Campbell and Haussdorff (BCH) [13). However, the exponential 
identities of BCH type involve infinite series and do not offer closed solutions. The key to the 
techniques of squeezing are the following formulae of Yuen [1| and Ma and Rhodes [2], which 
might be interpreted as exactly soluble cases of BCH and Zassenhauss. If no linear terms in H (t) 
are present, they read: 

U{t) = e B(tW\si{tWa e E(tW (Yuen, 1976) (7) 
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and 


u(t) = e l^(‘)« t2 -^(‘)« 2 l / 2 e -i«(tWa ( Ma and Rhodes , 1988) (8) 

where B(t), Q(t ), E(t ), z(t), «(() are c-number coefficients and = means the proporcionality of 
the unitary operators (U = U' => U = e in U’ } a € R). These identities precisely provide the proof 
that the variable oscillator potentials (more generally: quadratic, time dependent Hamiltonians) 
can produce the effects of squeezing (2) (or the scale transformation (3)). A particularly simple 
identity for two oscillator Hamiltonians Hj = p~ j 2 -\-u)jq~/2 acting during two time lapses Tj , t%, 
was detected by Griibl [ 1 4 J . If the time intervals Ti, To are in adequate proportion to the frequencies 
o,’i, U 2 (e.g., u)\T\ = 7r/2, U 2 T 2 = 37 t / 2) , then: 

e i\(qp+pq)/2 _ e -iT 2 H 2e -iTiHi , ( 9 ) 

where A = lna^/^i- (Note, that Griibl had no confidence to the operatorial formulae. He has 
proved (9) implicitely, working with Gaussian packets). 

While mathematically complete, (7-8) are not quite easy to apply, due to involved systems 
of non-linear equations for the c-number coefficients. This explaines a quick development of 
alternative methods derived from evolution matrices or adiabatic invariants [15, 16). Yet, the 
“damped oscillator” of 1940 [17], and the “step-Hamiltonian” of Griibl (9) remain the principal 
cases solved with all numerical details. 

It will be our purpose to show that the trend of the algebraic identities (7- 8) is not at all 
exhausted! To the contrary, appart of (7-8), it can provide a class of “easy formulae” for the 
squeezing and for more general control operations. 


2 The spin echo without spin 

The first “easy formula” (accidentally detected in 1977 (6[) has the form of the “circular identity”: 

*P 2 /2g— *9 2 /2g— tp 2 /2 > , e.~ iql l 2 = 1 (10) 

‘ 12 terms * 


for the operators q , p in L 2 (R), with [7, p) = i. Note, that all signs in the exponents are the same: 
the product (10) is therefore simpler than it could be in the classical case! The formula (10) has an 
elementary operational sense. Every operator expj — ip~/2\ represents the free evolution per unit 
time of the Schrodinger’s particle in L 2 (R). Every exp[ — iq 2 /2\ is the unitary evolution operation 
caused by infinitely sharp and quick pulse of time dependent oscillator potential (the “5-like kick” 


of the ellastic force): 


e^' 2 = lime- itLp2/2+< ‘ V/21 
€— ► 0 


(11) 
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The identity (10) describes a dynamical “evolution loop' 5 : the wave packet in L 2 ( R), manipulated 
by 6 oscillator kicks and 6 free evolution intervals must return to its initial state (no matter what 
this state was!). This might be illustrated by the following dosed diagramme: 


1 


I 



whose sides simbolize the free evolution intervals and vertices the oscillator kicks. An immediate 
consequence of (10) is: 


e +ip 2 / 2 _ e -iq 2 /2 e ~ip 2 /2 > # I e -ip 2 /2 e -iq 2 /2 

. 1 f i 


(12) 


The right hand side represesents a sequence of admissible dynamical events (6 kicks and 5 rest 
intervals), while the left one is the operator inverse to the free evolution. The formula (12) thus 
tells how to invert the free evolution. Since (12) is an operator identity, the prescription can be 
applied “in blind”: every wave packet in L 2 (R), entertained by 11 dynamical events must “go 
back in time”, returning to its past shape, no matter what this shape was [8], (Compare “Particle 
Memory” of Brewer and Hahn [18].) 

After some consideration, the formula (10) looses a part of mystery: it is just a “discrete 
imitation” of the oscillator force (the oscillator potential acts only in selected time moments, 
producing nonetheless a closed dynamical process). Note however the existence of other “circular 
identities” [8]: 

e -i P 2 / 2 e — , v / 3‘Z 2 /2 e — ip 2 /2 e +iv/3<7 2 /2 , . , e +i\/3g 2 /2 = (13) 

1 12 1 

The left hand side represents a sandwich of the 3 attractive and 3 repulsive pulses interrupted 
by 6 free evolution intervals. One might expect that the attractive and repulsive shocks will cancel 
“in average”, producing a zig-zag equivalent of the free evolution. However, it is not the case. The 
whole sequence traps the Schrodinger’s packets into a closed dance, with the evolution operator 
= 1. Note furthermore: 

(e -ip2 / 2 e -l<r ) = 1, (4 shocks, 4 free evolutions) (14) 
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Both formulae (13-14) can be illustrated by graphs: 



The simplest loop in L 2 ( R) must involve at least 3 kicks and 3 rest intervals; its general form is: 
e -np 2 /2 e — *«iV/2 e — »0 P 2 /2 e — *7r? 2 /2 e -m P 2 /2 e -^r q 2 /2 _ x (15) 



Its “incomplete version”: 

g +i7P 2 /2 _ e -iarq 2 /2 e -H3p 2 /2 e -iyrq 2 /2 e -iap' i /2 e -if)rq 2 /2 (16) 

permits one to manipulate the free evolution at will. Thus, for «,/?,7 > 0. (16) provides a 
prescription of how to enforce the Schrodinger’s wave packet to go back in time , wheieas for 
7 < 0, «, 0 > 0 one obtains a “time machine” able to slow or accelerate the free evolution [8, 10]. 

The loop formulae are the obvious analogue of the spin-echo for non-spin states. As far as 
we could check, the possibility of the (non-adiabatic) loop effects in L (R) was fiist piedicted 
in 1970 (by reinterpreting the transparency phenomenon of the potential wells; see Malkin and 
Man’ko [19], p. 388), though the subject was later pursued in a different direction. The first 
kicked system was considered in 1977 |6] and the manipulation of quantum states by potential 
pulses was systematically studied since 1986 |8, 9, 10]. 

The exponential identities suggest also how to generate the scale transformation. The simplest 
formula requires again a pair of oscillator pulses of different amplitudes: 


— tAp 2 /2g— i(l + l/A)g 2 /2g— i/> 2 /2g— i(l+A)</ 2 /2 = In A( 9 p+p 9 )/2 p ( 17 ) 

and produces the scale transformation superposed with parity (P). The repetition of the operator 
sequence of (17) yields the genuine squeezing (without parity: one of the simplest cases of Brown 
and Carson algorithm [20]). Some more general scenarios for the squeezing operation (2) (multiple 
kicks on a background of a constant ellastic force) are recently studied [21]. 
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3 Evolution control in three dimensions 


All these techniques concern the Schrodinger’s particle in 1 space dimension and are, in fact, 
only an abstract introduction to physically important problems. It is thus essential to find their 
analogues in 3 space dimensions. Some results can be already reported. 

In the first place, the sequences of sign changing kicks [e.g. (13-14)) can be used to construct 
sequences of harmonic pulses in R 3 generating the loop effect in L 2 (R 3 ) [8], This suggests, that 
the loop effect (state echo) in R 3 might be produced, in principle, by shock waves of source free 
external fields. As the matter of fact, some closed dynamical processes can be induced even without 
any kicks, by a source free, stationary field of an adequately gauged ion trap |9). A simple scenario 
of additional potential kicks (electric pulses applied to the trap walls) permits then to generate 
effects of squeezing upon the charged wave packet retained in the trap interior (see the report by 
one of us |9]). 

What no less important, the effects of positive (attractive) oscillator potentials in L 2 ( R) tra- 
duce themselves immediately into effects of homogeneous magnetic fields in 3 space dimensions. 
As an example, we have considered a quite simple sequence of identically shaped magnetic pulses 
in three orthogonal directions: n, m, s, n, m, s • • •. As we have reported on the previous IWSSUR 
93, an adequate proportion between time separations and the pulse intensity assures that the 
sequence must produce the loop effect for the wave packets in L 2 (R 3 ). Moreover, the same opera- 
tional scheme, with differently shaped pidses, turns out to work as a “time machine” , permitting 
to accelerate, slow or invert the free evolution operation of the Schrodinger’s wave packet [10]. 



A sequence of homogeneous magnetic pulses from 3 orthog- 
onal directions permits to manipulate the free evolution (see 
our report in IWSSUR 93), 

4 The general manipulation scheme 

The most immediate reason why one might be interested in the “evolution loops” is the possibility 
of controlling the fuzziness (diffraction of the wave packets due to its free evolution), essential for 
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electronic microscopy, programming the non-demolishing measurements etc. (see also C aves et 
al. [22], Yuen [23|, Royer [24]). Also, the original subject of the transparent wells [19] might 
still bring some surprises |25|. However, the loop phenomenon seems most crucial for the general 
manipulation methodology. 

The class of the dynamical operations induced by stationary fields is rather narrow (foi the 
Schrodinger’s particle they are always of the form exp[-i//| where the exponent H is at most 
quadratic in p!). 

The situation is more interesting for a microobject trapped in an oscillating field of an evolution 
loop. As long as the loop fields are mantained, the wave packets perform a “periodic dance” . A 
distinct phenomenon occurs, if the loop fields are perturbed or imperfect. Instead of a closed 
process, the system will then perform, after every loop period r, a non trivial unitary operation, 
interpretable as the loop precession. 



The precession of a distorted loop: 
a natural key to the manipulation. 

An elementary algebraic argument shows that the precession operations are much more general 
than the operations stimulated by the stationary fields. In fact, they are the key to solve the 
manipulation problem: by “adding precessions” an arbitrary unitaiy opeiation U . Ti, ► be 

approximated [8j. In some cases, already an unsophisticated distortion of the circular processes 
brings interesting results (like e.g. the squeezing or free evolution distortion in “wrong loops”). 
In principle, every one of the “circular identities” (10,13-15) is a natural staiting point foi some 
manipulation procedures solving the inverse evolution problem (1). With one little amrnendment, 
however. 

5 The “soft kicks” 

The “ellastic kicks”, while of undeniable illustrative value, are not so easily accessible in labora- 
tories. The difficulty is almost anecdotic if the “<5-like kick” has to be engineered with the help of 
homogeneous magnetic field acting e.g. inside of a cillindric.al solenoid. Since of the lesulting 

“magnetic oscillator” is proportional to B(t) 2 [10], (where B(t) is the magnetic field intensity), 
the B (t) in the solenoid would have to model the square root of the Dirac s 6. The lequest might 
be promising for the theory of non-linear distributions, but is a nightmare in the laboratory! (As- 
suming even that the laboratory team would dominate the techniques of approaching y/sit), the 


247 


radiative corrections would probably spoil the effects of the operation). 

What one needs are the soft analogues of oscillator kicks (11), and they are not so difficult 
to programme with the help of exponential formulae. Below, we shall report a quite simple 
“exponential experiment” . 

Consider first of all the product of three operators: 


W = e iX ( qp+pq ^ 2 e~ ilql I 2 e _ia (P 2+ 9 2 )/ 2 


( 18 ) 


Let’s ask the question: can one choose A, 7,0 to be three functions of time in such a way that 
the product W fulfills a physically interpretable evolution equation 


dW 

~IT 




(19) 


with H ( t ) having the oscillator form ( 6 )? To simplify the problem, we shall first determine A and 
7 as functions of a, A = A(o) and 7 = 7 ( 0 ), and only afterwards we shall look for a = a(t). Each 
term in (18) is easily differentiable: 


dW 

da 


HoMc- 1 ^) W („) 


where H 0 = p 2 /2 + q * / 2 = a^a T 1/2. Due to the transformation rule (4) and: 




(? + 7 <?) 2 r 

2 2 


one easily finds: 

dW 


* — = H(a) W(a) = ( (-A + 7)^^^ + + (7 + 7 2 + \)^ q ~ ) W(a) 


da 


2 


To assure that the term with (qp -f pq)/‘2 vanish it suffices to put: 

7(0) = A(«) 

thus obtaining: 

H{a) = e~ 2xP j + (7 + 7 2 + l)e 2A y 
If now a = a(t), the differential equation for W in terms of t reads: 


dW (da\~ , 

ir “hr 


( 20 ) 


( 21 ) 

(22) 

(23) 

(24) 

(25) 
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and if in addition o(f) is determined by: 


da _ c ~‘2Mn) 

(26) 

(it. 


then (25) acquires the familiar form: 



(27) 

where a(£) is given by (26) and 


g{t) = (7 + T + l)c 4A(ft(,)) 

(28) 


If now W(t) satisfies (27), then U (t) = W{t)W(0) 1 solves the evolution problem ( 6 ) with 
< jj[t ) 2 = g(t) and with the initial condition f 7 ( 0 ) = 1 . 

Adopting A(o) defined in [0, 27r] as our arbitrary “manipulation function”, we can model at will 
the desired properties of W (t) [and consistently, of Thus, if A(0) ^ A(27r) but A(0) — A(27r), 

(/ at a = 27 r becomes the scale operator of form (4) and the function g(t) defined by (28) gives 
the prescription of how the effect can be generated. If, however, A(0) = A(2 tt) but 71 = A(2tt) / 
A(0) = 0, the same product (18) reduces to the single non trivial term 

W(t) = e" i7,,a/2 (29) 

and henceforth, U imitates the effects of the 5-like kick of the oscillator force. As an example, we 
report two simple computer simulations where the manipulation function A(o) yields either the 
“soft imitation” of the oscillator kick or the coordinate squeezing (see below). It seems pertinent to 
notice, that if an authentic kick were to be applied in the laboratory e.g. by creating a very short 
and sharp magnetic pulse, then in the first place it could never be exact (nor well approached: 
the 5 -functions are not truly accessible in labs!) In our scenario below, this difficulty is absent: 

Manipulation function Ellastic amplitude g(t) 


A.(a) g(t) 



^(a) = sin(a/2) - sin(a)/2 Effective Operation: 

e“ |2/2 (repulsive kick) 
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The kick effect can be exact (produced with unlimited accuracy), even if negative, and is 
achieved by softly varying fields awaking little radiative response. By the same, all previous results 
involving the oscillator kicks [6-9,23] can be interpreted as realistic laboratory prescriptions. Note 
also the squeezing scenario based on the same formula (18): 


Manipulation function 

k(a) 



Ellastic amplitude g(t) 

9(t) 

2 

1 

0 2 4 6 



a.(a)=0.05 cosa-02 cos(a/2) Effective Operation: 

e K0.4)(pq-qp)/2 

The shape of g(t) agrees with the observation that the squeezing is caused by an increase of the 
ellastic constant [14, 17]. 

The story does not end up here; it hardly starts. The method of distorted loops makes 
possible much more sophisticated manipulations of quantum degrees, which will be probably the 
daily routine of the experimental physics in a predictable future. 
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